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A  FORTRAN  IV  PLOT  ROUTINE  WITH  HIDDEN 
LINE  SUPPRESSION  FOR  USE  WITH  SMALL  COMPUTER  SYSTEMS 


INTRODUCTION 

A  frequently  encountered  problem  in  computer  graphics  is  the  presenta¬ 
tion  of  plots  representing  the  value  of  a  single  dependent  variable  as  a 
function  of  two  independent  variables.  A  commonly  used  solution  to  this 
problem  consists  of  a  depiction  of  a  grid  located  on  the  surface  represent¬ 
ing  the  function.  The  curves  comprising  this  grid  are  obtained  by  letting 
the  remaining  independent  variable  vary.  The  viewpoint  is  then  established 
by  coordinate  rotation,  with  respect  to  the  plane  of  the  graphics  medium, 
and  the  projection  of  the  curves  comprising  the  grid  are  plotted  on  this 
plane.  In  order  to  avoid  a  confused  or  ambiguous  plot,  some  method  must  be 
used  to  suppress  portions  of  the  curves  which  would  be  hidden  from  the  view¬ 
ers  position. 

This  report  presents  a  routine  which  has  been  found  useful  in  solving 
this  hidden  line  problem.  Simplification  of  the  general  problem  has  been 
used  to  provide  a  method  which  is  suitable  for  use  with  small  computer 
systems.  The  routine  was  developed  for  use  with  a  Data  General  Corporation 
NOVA  800  with  associated  Tektronix  Incorporated  graphics  terminal  equipment. 
The  coding,  however,  is  Fortran  IV  and  is  easily  adapted  for  use  with  other 
computer  systems. 

DESCRIPTION  OF  THE  ROUTINE 

Plots  made  using  the  routine  are  shown  in  figures  1  through  4.  A  list¬ 
ing  of  the  Fortran  IV  source  code  is  included  as  Appendix  A.  The  comments 
included  in  the  source  listing  should  allow  the  program  flow  to  be  followed 
once  the  parameters  and  general  method  are  explained. 

Consider  the  function  to  be  represented  to  be  of  the  form:  y  -  F(x,z). 
We  wish  to  depict  curves  of  y  ■  F(x  •  const.,  z)  and  y  *  F  (x,z  •  const.). 
Taking  one  of  these  curves,  the  coordinates  are  transformed  by  a  rotation  of 
<1>  in  the  horizontal  plane  (x,z  plane)  and  then  tipped  forward  by  an  angle,  9. 
The  projection  in  the  plot  coordinates  is  then  obtained.  The  angle  $  is 
required  to  be  45  degrees.  This  allows  the  entire  problem  to  be  handled  in 
terms  of  indices  of  arrays,  which  simplifies  the  program  and  economizes  on 
storage  requirements.  The  coordinate  transformations  are  the  responsibility 
of  the  calling  program.  Listings  of  two  routines,  useful  in  this  regard, 
are  included  in  Appendix  B. 

The  transformed  data  is  stored  in  two  arrays  (IXPA  and  IYPA)  giving  the 
projected  coordinates  of  the  data.  Hidden  line  suppression  is  accomplished 
by  comparison  with  an  array  (IYMA)  which  stores  the  maximum  Y  encountered 
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as  a  function  of  index.  IYMA  must  be  properly  initialized  in  the  calling 
program  but  is  subsequently  maintained  by  the  plot  routine.  Plots  should 
be  accomplished  from  front  to  back.  That  is,  constant  z  cuts  must  start 
with  maximum  z  while  constant  x  cuts  must  start  with  maximum  x.  Furthermore, 
at  the  completion  of  one  set  of  variable  equals  a  constant  cut,  IYMA.  must  be 
properly  reinitialized  to  the  projected  view  of  the  front  of  the  surface 
before  proceeding  to  plot  the  other  set. 

Provision  is  made  for  updating  of  the  IYMA  array  for  values  of  the 
independent  variable  between  plotted  constant  cuts.  This  gives  a  truer 
picture  of  the  surface  as  a  two  dimensional  solid  (see  figures  1  and  3)  but 
may  not  be  desired  in  some  Instances.  For  example,  if  either  constant  x  or 
constant  z  cuts  are  used  alone,  the  presentation  obtained  by  not  updating  for 
intermediate  surface  may  be  preferred  [1]  (also  see  figure  4).  For  experi¬ 
mental  data,  intermediate  values  for  the  "surface"  are  frequently  unknown 
and  consequently  demand  the  use  of  this  option.  The  routine  calculates  the 
intersection  of  the  line  segment  joining  adjacent  data  points  and  the  line 
segment  joining  the  pertinent  stored  maximum  y  points  and  only  plots  the 
appropriate  section  of  the  functional  curves.  Two  plot  routines,  MOVE  and 
DRAW,  are  referenced  in  the  routine.  Since  these  are  hardware  dependent, 
listings  are  not  Included,  but  functional  equivalents  must  be  supplied  by 
the  user.  MOVE  (IX,  IY)  is  equivalent  to  a  "pen  up"  command  while  DRAW  (IX, 
IY)  is  equivalent  to  "pen  down",  both  with  plotter  motion  to  the  point  (IX, 
IY). 

DESCRIPTION  OF  PARAMETERS 

The  call  to  the  routine  is: 

CALL  P3D  (IXPA,  IYPA,  IYMA,  MDIM,  NDIM,  NPTS,  MLOW,  IPL) 

Arrays  are  dimensioned: 

IXPA  (0:NDIM)  ,  IYPA  (0:NDIM),  IYMA  (MDIM: NDIM) 

Definitions  of  parameters: 

IXPA  is  the  x  plot  array. 

IYPA  is  the  y  plot  array. 

IYMA  is  the  array  storing  maximum  value  of  y  encountered  for  a 

particular  index. 

MDIM  should  be  set  "  -NDIM  in  the  calling  program 

MLOW  must  be  set  equal  to  the  lowest  value  of  index  for  IYMA 

(i.e.,  the  left  most  point)  in  the  calling  program. 

NPTS  is  number  of  points  to  be  plotted. 

IPL  *  0  don't  plot,  but  update  IYMA 

1  plot  and  update  IYMA 

The  relationship  of  the  various  indices  is  diagrammed  in  figure  5.  This 
figure  is  a  view  looking  down  on  the  original  coordinate  system.  Note  that 
$*45°,  variation  in  9  does  not  effect  these  relationships.  This  is  the  key 
to  the  simplicity  of  the  hidden  line  routine. 


2 


SUMMARY 


A  routine  useful  for  plotting  "three-dimensional"  representations  of 
single  valued  functions  of  two  independent  variables,  or  the  equivalent  in 
experimentally  obtained  data,  has  been  described.  The  method  used  limits 
the  rotation  of  point  of  view  in  order  to  economize  on  required  computer 
resources,  thus  making  the  method  applicable  to  small  computer  systems. 
Hidden  line  suppression  is  implemented,  including  the  optional  suppression 
of  those  portions  of  plotted  cuts  which  would  be  hidden  by  the  surface 
between  cuts,  in  the  event  these  intermediate  values  for  the  functions  are 
known. 

REFERENCE 

[1]  "A  CDC  3800  Computer  Subroutine  for  Producing  Isometric  Plots  on  the 
565  Cal  Comp  Plotter",  H.A.  Brown,  NRL  Computer  Note  31,  1  Aug  1969. 


Fig.  1  -  Plot  of  a  surface  using  both  constant  x  and  constant  z  cuts. 
The  YMAX  array  is  updated  between  plotting  cuts  to  obtain  an  accurate 
representation  of  the  surface. 
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Fig.  2  -  Plot  of  surface,  as  in  figure  1,  except  the  YMAX  array  is  not 
updated  between  the  plotted  cuts.  The  deficiencies  in  the  representation 
of  the  surface  are  apparent. 


-  Plot  of  surface  using  constant  x  cuts  alone.  The  YMAX 
array  is  updated  between  plotted  cuts. 
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Fig.  5  -  Diagram  showing  Che  relationship  between  Indices  of  the 

various  arrays 
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Appendix  A  —  Listing  of  the  plot  routine 


TVf*  ^SUIROUTINE  P3D( IXPA, IYPAt I VNA.HDIR.HDIN, NPTS.  fJLOU.  IPl ) 

C  ROUTINE  TO  PLOT  FROR  ARRAYS  SIANKIM  HIDDEN  LINES 

C  NOOIFIE*/ TO  INTERPOLATE  PARTIALLY  HIDDEN  VECTORS 

C  upC  9/1/11 

DIRENSION  IXPA(#lNDIH  >.IYPA(#IHDIN), IYNA<HDIN«NDIN) 
IFUPL.CO.O)  QO  TO  41  jftOH'T  PLOT.  BUT  UPDATE  RAX  ARRAY 
DO  31  IN0«1.NPTS 
I-1ND-1 

NaRLOU+I{ INDEX  OF  RAX  ARRAY 
IF(I.NE.S)  00  TO  33 

CALL  HOUE(IXPA(S),IYPA(S) IjPOINT  IS  lEOINNINQ  POINT 
GO  TO  31 

33  IVl-IVPA(I-l) 

IY8-IYPACI) 

IYNWYHACR-l) 

ivn8«!vnA(tn 

XXlaIXPA(l-l> 

ixa*ixPA(i> 

C  ARE  BOTH  PTS  ABOVE  (OR  BE LOU)  YRAXT 

IF( .NOT. ( IY1 .LE< IVR1 .AND. IY3.LE. IVH2 ) >  00  TO  SI 
CALL  H0UE(IX2.IY2)  | BOTH  PTS  BELOU 

GO  TO  31 

St  IFf .NOT. ( IV1.GT. IVN1 .AND.IY2.CT.IVN2 ) )  00  TO  S2 

CALL  DRAU( 1X2, IY2)  jBOTH  PTS  ABOVE 
GO  TO  31 

C  ONE  END  OF  VECTOR  IS  HIDDEN.  SET  UP  TO  INTERPOLATE 

52  sn«FL0AT(IYn2-IYHl)/(IX2-IXl)  .SLOPE  OF  YNAX 

SP*PLOAT( IY2-IY1 )/(IX2-IXl )  , SLOPE  OF  PLOT 

X*IXl«fLOAT(IYl-Ivni )/(SN-SP)  .X  AT  PT  OF  INTERSECTION 

v-sp»(x-float<  ixi  mm  .v  at  pt  of  intersection 

IX*IFIX(X)  | INTEGER  COORDINATES 
IY.IFIX(Y)  t0F  INTERSECTION 
C  OCTERNINE  WHICH  END  OF  VECTOR  IS  HIDDEN  AND  PLOT 

IF( .NOT • (IV1 .GT. 1YN1 -AND. IV3.LE. IYM2 ) )  GO  TO  S3 
CALL  DRAU(IX.IY)  jLEFT  ABOVE 
CALL  N0VE(IX2. IY2)jAND  RIGHT  BELOU 
GO  TO  31 

53  CALL  HOVEdX.IY)  jLEFT  BELOU 

CALL  DRAUdX2,IY8)|AND  RIGHT  ABOVE 

31  CONTINUE 

41  DO  35  INDal,NPTSj UPDATE  NAX  ARRAY 

I-IND-1 
NaNLOU+I 

3S  IF< IYPAC l ).GT . lYHA(R))  1YNA<N  )-IYPA( I ) 

RETURN 

END 
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Appendix  C  —  Listing  of  program  used  to  generate  plots  shown  in  the  figures 

C  DOM  PROGRAM  FOR  HIDDEN  IXNC  PLOTTING,  HPC  D/S/D  1 

dimension  iypa<*iim>.ixpa<*sisg).ivha<-i6*iis*) 

CALL  XNITT 
IS-* 

PI-3.14159B? 

R  AD-1 1*. /PI 

ACCEPT*  IXC7  *.IXC  ) CENTER  OP  SCREEN  COORDINATES 


99 


ACCEPT*  I VC?  *»IYC  i  • 

ACCEPT*  NMAX7  *.NMAX  .MAX  LENGTH  OP  AXES  IN  SCREEN  POINTS 
PNNAX-NNAX  > FLOAT  AXIS  LENGTH 
NOIV-IS*  .NUMBER  OF  POSS1DLE  DIVISION  ALONG  AXES 
NPTS-NDlU+1  .NUMBER  OF  PLOT  POINTS 

ACCEPT*  NCUTr  *,NCUT  .NUMBER  OF  DIVISIONS  ON  AXES  DETUEEN  PLOTTED  GRIDS 

ACCEPT*  UPDATE  NAX  ARRAY  DETUEEN  GRIDS')  (TYPE  1  FOR  YES»  *,NS 

NDIN—  IS*  .SET  LOWER  DINENSION  OF  VNAX  ARRAY 

NDIN-1S*  ;SET  UPPER  DIPENSIONS  OF  ARRAYS 

TYPE 

TYPE*  AT  BREAKS  WITH  SOUND  OF  DELL* 

TYPE*  RESUNE  OPERATION  IV  INPUTTING  NUNIER  AND  CR* 

TYPE*  (INPUT  OF  1  RESTARTS  PROGRAM  AT  BEGINNING)* 

TYPE 

ACCEPT*  INPUT  1  FOR  PLOTTER  OUTPUT  *,NP  j SELECT  SCREEN  OR  PLOTTER 
IF(NP.NE.l)  GO  TO  SS 

TYPE*  SUITCH  TO  AUX  OUTPUT,  MCE  TO  TTY  AT  END  OF  PLOT* 

PAUSE 
CALL  PLON 

IF(HP.NE.ll  CALL  CLEAR  ; CLEAR  SCREEN  FOR  GRAPHICS 


CO  TO  PLOTS 


29 


NLIN-2*NDIV*l 
DO  29  N-l.NLIN 
I-N-NDIV-1 
IYNA(I )■* 

CALL  GRAPH(9€) 

DO  Z-CONST  CUTS 

DO  37  N-l.NPTS 
J-N-l 


I CLEAR  NAX  ARRAY 
I SET  GRAPH  NODE  AND  VECTOR  TYPE 

I LOOP  OVER  Z 


Z-FNNAX-JKFNNAX/ND1V)  «Z  COORDINATE  OF  CUT 
PLOU-PDIP-J  |SET  LEFT  POST  POINT  OF  NAX  ARRAY 
IFjnODU.NCUTi.EO.O)  GO  TO  34  , DETERMINE  VHETHER  TO  PLOT 
IPL-* 

GO  TO  36 
34  IPL-1 

^36  IF(NS.EO.*.AND.IPL.EO.*>  GO  TO  37  |DO  NOT  UPDATE  FOR  INTERMEDIATES 

0°  33  N-l.NPTS  | PUT  TRANSFORMED  X,Y  IN  ARRAYS  AM)  PLOT 
I-N-l 

X-It(FNNAX/MDIV> 

TX-X-FNNAX/3.  > TRANSLATE 
TZ-Z-FNNAX/R. 

RHO- SORT  t  TXf t2-TZtta ) 

ARO-(RHO/(FNNAX/|.))«S.tPI/«. 

V-FNNAXaADSlOOSl ARG )  UCOSIRHODPI/FNNAX ) 

IFCRNO.OT, (FNMAX/I, ) )  V«9. 
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Appendix  C  (Cont’d)  —  Listing  of  program  used  to  generate  plots  shown  in  the  figures 


33 

37 


.trams mr 

cau  p39<ixpa7xypa.iyha.roin1mir<npts,niou,xpl>  »plot  array 

CONTINUE 


CALL  HQUCCO.lO)  (ROUE  HEM  DOTTON  LEFT 
CALL  ALPHA  (3)  i RETURN  TO  ALPHA  (WOE 
IF(HP.EQ.l)  CALL  PLOFF 


CALL  KLL 
ACCEPT  ICON 

IFCXCON.EG.l)  00  TO  99  |  RESTART 

XF(NP.EO.l)  CAU  PLON 

CALL  GRAPH <  96)  j  RESET  GRAPH  NODE  AW  UCCTOR  TYPE 


oo  x-coMsr  curs 
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00  39  N-l.NUN 
I-N-NDIU-l 
IVHAI X  >*• 


I  RECLEM  RAX  ARRAY 


l LOOP  OVER  X 


46 


00  47  fl-l.NPTS 
J*H-1 

X«FNNAX- JS < FNNAX/MOIU )  »X  COORDINATE  OF  CUT 
NLOU— J  i  LEFTMOST  POINT  OF  IMX  ARRAY 
IF<NOD(J.NCUT)  EG. 9)  GO  TO  44  (PLOT  OR  UPDATE? 
IPL-G 
GO  TO  46 
IPL-1 

IFCNS.EQ.O.AND. 


_ XPL.EQ.S)  GO  TO  4?  iDO  NOT  UPDATE  FOR  INTERNEDXATES 

00  43  N-l.HPTS  (PUT  TRANSFORRED  COORDINATES  IN  ARRAYS  AMD  PLOT 
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I-N-l 

Z'FNNAX-IXFNHAX/MXU) 

TX*X-FNNAX/R . 

TZ*Z-FHRAX/8. 

RHO*SGRT(TX»|«TZ«ka  > 

ARC*CRHO/(FNAAX/a. USS.SPI/R. 

Y-FWIAXSASS  ( COS  ( AR6  )  >SCOS  ( RHOXPI^FNRAX ) 
IFCRHO.GT. (FNNAX/3. ))  Y»0. 

CALL  ISOPRIX.Y.Z.XP.YP.IS) 

IXPA! X )*XP4IXC 
IYPAUJ-YPaIYC 

CALL  P3DIIXPA.IYPA,IYNA,NDIN,N0IR,NPTS,nL0U,IPl> 
CONTINUE 

CALL  NOUE(O.O)  jNOUE  TO  LOWER  LEFT 
IF(NP.EG.l)  CALL  PLOFF 
CALL  ALPHA! 3 I  i RETURN  TO  ALPHA  NOOE 
CALL  SELL  (SIGNAL  EW  OF  PLOT 
ACCEPT  ICON 

IFdCON.EO.I)  GO  TO  99  (RESTMT 
*•♦0 
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